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LONG-TERM  GOALS 

My  long-term  goal  is  to  contribute  to  understanding  of  the  nonlinear  dynamics  of  the  wave  sea  surface 
excited  by  wind.  Of  particular  interest  to  me  is  the  role  of  such  interaction  in  long-term  prediction  of 
wave  amplitudes  excited  in  rough  seas. 

OBJECTIVES 

1  wish  to  develop  the  phenomenological  diffusion  model  of  interaction  of  gravitational  waves  on  water 
surface  in  presence  of  wind  and  viscosity.  Motivation  for  development  of  such  a  model  is  the  fact  that 
numerical  solvers  based  on  Hasselmann  kinetic  equation  for  waves  are  time-consuming  and  hardly  can 
be  used  for  practical  purposes.  Numerical  solver  based  on  the  diffusion  model  is  expected  to  be  at  least 
three  orders  of  magnitude  faster. 

APPROACH 

1  propose  two  generations  of  simplified  model  of  Hasselmann  kinetic  equation  for  waves.  The  first 
generation  model 
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Eq.(l)  is  the  diffusion  approximation  obtained  from  differential  approximation  to  the  kinetic  equation 
[2],  [3].  One  should  emphasize  that  equation  (1)  is  the  model  equation  and  the  value  of  the  constant 
a  can  be  defined  only  from  comparison  with  the  results  of  numerical  simulation  of  the  kinetic  equation 
or  laboratory  measurements.  The  results  of  numerical  simulation  of  the  Eq.(l)  have  shown  good 
agreement  of  angle  averaged  frequency  directional  distribution  of  the  nonlinear  interaction  term  with 
the  corresponding  results  for  kinetic  equation  (see  [4],  [5]).  The  correspondence  of  angular 
distributions  for  particular  frequency  values  was  not,  however,  very  good. 
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The  second  generation  diffusion  model  takes  into  account  the  effects  of  non-locality  in  the  nonlinear 
interaction  term  while  preserving  such  important  properties  of  the  first  generation  diffusion  model  as 
scaling,  conservation  laws  and  Kolmogorov  solution: 
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and  {n^'^  = 

Besides  the  term  in  the  right-hand  side  that  is  similar  to  the  first  generation  model  there  are  the 
other  terms  representing  all  possible  combination  of  the  terms  which  are  cubical  with  respect  to  n  and 
contain  an  integration  over0  .  We  intentionally  omitted  the  term  (j  ndcp)  and  I  n^dcp  as  they  make  no 
contribution  into  angular  anisotropy.  As  in  the  case  of  first  generation  diffusion  model,  the  coefficients 
=  1,4  are  defined  Ifom  the  comparison  with  the  results  of  numerical  simulation  of  kinetic  equation 

or  experimental  measurements.  Obviously,  the  choice  of  the  coefficients  Wj  ^  0,a2  =  =0,a^  =0 

corresponds  to  the  first  generation  model. 


where  (^n)  =  —  j  ndcodcp 


The  following  is  the  description  of  the  algorithm  of  definition  of  the  set  of  coefficients  a^,i  =  1,4  .  Let 
S„,{n{co,(p))  be  nonlinear  interaction  term  at  the  kinetic  equation,  rij  =  nj{co,(p) 

4 

be  JONSWAP  initial  conditions  [4],  [5]  and  S„,  =  (n(ft),^))be  nonlinear  interaction  term  at  the 

i=i 

second  generation  model  (3),  where 

5'j(n)  =  Ln^  ,S2in)  =  Ln^  in'),S2{n)  =  Lnin^'^,S^{n)  =  Lninf 


We  are  looking  for  the  best  approximation  to  the  function  S„i{nj)'by  the  linear  combination  of  the 
functions  S^  {rij  ),  i  =  1,4  in  the  norm  : 


da, 
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The  coefficients  a,  ,  /  =  1,4  are  the  solution  of  the  system  of  equations  (4)  which  is  the  system  of  linear 
equations 


ByCCy  =  Cj,i,j  =  1,4 

where  symmetric  matrix  and  vector  Cy  are  defined  by 

By  =  \  SXnj)S  j{nj)dcod(p,C  j  =  \  S„i{nj)Sj{n  j)dcod(p,i,  j  =  \,A 


(5) 


Both  the  first  and  second  generation  diffusion  models  overestimated  the  nonlinear  interaction  in  the 
case  of  narrow  directional  spectra.  In  FY2000,  to  fix  this  a  third  generation  model  was  formulated  in 
which  there  is  an  expansion  about  the  wind  direction.  For  very  narrow  angular  spectra,  this 
formulation  should  vanish  as  expected  Ifom  the  full  calculation.  Scaling  coefficients  must  be  found  by 
comparison  t  the  exact  Hasselmann  equation. 

Additionally,  we  proposed  a  “clean  test”  for  checking  the  accuracy  of  the  various  nonlinear 
approximation  models.  This  calculation  is  based  on  control  of  the  conservation  constants;  action, 
momentum  and  energy. 

In  a  spearate  action,  a  model  for  rogue  waves  was  also  proposed  based  on  the  integral  equation 
containing  the  ouplinc  coefficient  for  surface  waves.  By  using  a  simple  approximation  for  the 
coefficient  based  on  the  modulus  of  the  interacting  wave  numbers  it  is  possible  to  arrive  a  a  quasi- 
soliton  solution  for  the  formation  of  localized  wave  groups.  Numerical  simulations  in  one  dimension 
suggest  that  quasi-solitons  can  merge  in  non-elastic  collisions  generating  waves  of  very  high  amplitude 
from  initially  smooth  conditions. 

The  numerical  part  of  the  work  was  performed  by  Dr.  Andrei  Pouchkarev,  Department  of 
Mathematics,  University  of  Arizona. 

WORK  COMPLETED 

First  and  second  generation  diffusion  model  were  proposed.  Based  on  these  models  supplied  by 
viscous  and  forcing  terms,  computer  codes  have  been  created.  Both  codes  have  been  shown  to 
reproduce  important  properties  of  the  original  Hasselmann  equation  —  stationary  Kolmogorov 
solutions.  Typical  calculation  of  the  evolution  of  the  turbulence  driven  by  external  forcing  from  the 
initial  low-level  random  noise  conditions  to  the  stationary  equilibrium  state  on  the  grid  of  40x130 
nodes  in  angle-frequency  space  took  the  time  of  the  order  of  dozens  of  minutes  using  Pentium- 133 
MHz  PC.  Comparison  of  the  results  of  exact  Hasselmann  equation,  first  and  second  generation 
diffusion  models  has  been  made.  It  was  shown  what  are  the  advantages  of  the  second  generation  model 
with  respect  to  the  first  generation  one. 

RESULTS 

We  calculated  the  coefficients  a,  ,  /  =  1,4  as  the  solution  of  the  system  of  linear  equations  (5)  using  the 
functions  rij  and  from  numerical  simulation  of  kinetic  equation  for  waves  by  Resio  and  Tracy 

[4],  [5]:  =-2.35-10 ',a2  =5.14-10  ",a3  =-9.76-10  ",a4  =4.67  10  ‘. 

Both  first  and  second  generation  model  qualitatively  correctly  reflect  the  behavior  of  (/,^)in  the 
kinetic  equation.  It  should  be  noted  that  both  of  them  fail  to  reproduce  such  fine  details  of  the 
Sni{f,(p)  distribution  as  two-hump  maximum  and  small-amplitude  high-frequency  tail. 

Fig.  1  represents  directional  (angle  averaged)  distribution  of  the  spectrum  as  a  function  of  the  frequency 
for  three  models.  The  correspondence  of  the  second  generation  model  with  Resio-Tracy  model  is  better 
than  first  generation  model. 


Fig. 2  represents  angular  dependenee  at  the  frequeney  of  the  minimum  of  the  speetrum  $f=0.16$. 
Seeond  generation  model  improves  eorrespondenee  with  Resio  -  Traey  model  with  respeet  to  the  first 
generation  model. 

Fig.3  represents  nonlinear  interaetion  term  distribution  as  the  funetion  of  the  Ifequeney  at  zero  angle. 
While  the  seeond  generation  model  dependenee  is  mueh  eloser  to  Resio-Traey  results  in  the  vieinity  of 
the  minimum  and  the  maximum  of  the  speetrum,  it  still  fails  to  describe  high  -  Ifequeney  tail  of  the 
distribution. 


Figure  1.  Directional  (angle  averaged)  spectrum  as  a  function  of  the  frequency  f . 

Crosses  represent  Resio-Traey  results,  dashed  line  —  first  generation  model,  solid  line  - 

second  generation  model. 


Figure  2.  Angular  dependence  at  fixed  frequency  f  =  0.16.  Crosses  represent  Resio  - 
Tracy  results,  dashed  line  -  first  generation  model,  solid  line  -  second  generation  model 


Figure  3.  Nonlinear  interaction  term  dependence  on  frequency  at  the  angle  (p  =  0. 
Crosses  represent  Resio-Tracy  results,  dashed  line  —  first  generation  model,  solid 

line  -  second  generation  model. 


IMPACT/APPLICATION 

The  second  generation  diffusion  model  is  the  simple  model  which  takes  into  account  the  effects  of 
non-locality  of  the  interaction  of  surface  gravitational  waves.  As  the  first  generation  model,  it 
preserves  the  constants  of  the  motion,  has  righteous  scaling  and  Kolmogorov  solutions.  Second 
generation  model  also  contain  first  generation  model  as  a  special  set  of  parameters. 

Least  square  optimization  allows  to  find  unknown  coelficients  for  the  second  generation  model  using 
the  results  of  kinetic  equation  numerics.  Comparison  of  the  second  generation  model  with  kinetic 
equation  numerics  shows  that  second  generation  model  improves  directional  and  angular  dependence 
of  the  spectrum  with  respect  to  first  generation  model.  It  still  fails  to  describe  fine  details  of  the 
spectrum  such  as  two  -  hump  shape  of  the  maximum  and  the  high-frequency  tail. 

TRANSITIONS 

Our  first  and  second  generation  codes  have  been  handed  to  people  at  the  Army  Water  Experimental 
Station  for  testing. 

RELATED  PROJECTS 

There  are  two  related  projects: 

1.  "The  response  of  wind  ripples  to  long-surface  waves  —  application  to  radar  studies",  ONR  N00014- 
98-1-0439. 


2. "Statistical  model  for  deep  and  shallow  water  waves",  ARO  DACA  39-99-C-0018. 
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